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ABSTRACT 

Aims. The abundances of Fe in the intracluster medium of nearby (2 < 0.08) clusters were measured up to 0.3~ O.Sriso- 
Methods. We analyzed 28 clusters of galaxies observed with XMM-Newton. We derived Fe abundances from the flux 
ratios of Fe lines to the continuum within an energy range of 3.5-6 keV to minimize and evaluate systematic uncer- 
tainties. 

Results. The radial profiles of the Fe abundances of relaxed clusters with a cD galaxy at their X-ray peak have similar 
slopes. These clusters show similar enhancements in the Fe abundance within 0.1ri8o, and at 0.1-0.3ri8o, they have 
flatter Fe abundance profiles at 0.4~0.5 solar, with a small scatter. Most other clusters, including merging clusters, also 
have similar Fe abundance profiles beyond O.lrigo. These clusters may have universal metal enrichment histories, and 
a significant amount of Fe was synthesized at a very early stage in cluster formation. Mergers of clusters can destroy 
the central Fe peak. 
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1. Introduction 

Clusters of galaxies are the largest gravitationally bound 
objects in the universe. The intracluster medium (ICM) 
contains a large amount of metals, which are synthesized 
mainly by supernovae (SNe) in early- type galaxies (e.g., 
Arnaud et al. 1992; Renzini et al. 1993). Thus, metal abun- 
dances in the ICM provide important clues for understand- 
ing the metal-enrichment history and evolution of galaxies 
in clusters. Because both SN II and SN la synthesize Fe, the 
distribution of Fe in the ICM contains information about 
the star-formation histories of massive stars and the history 
of chemical-enrichment attributed to SN la. 

The ASCA satehite (Tanaka et al. 1994) first enabled 
us to measure the distribution of Fe in the ICM (e.g., 
Fukazawa et al. 2000; Finoguenov et al. 2000; 2001). The 
Fe abundances of these clusters are 0.2-0.3 solar, adopt- 
ing the solar abundance from the "photospheric" values 
given by Anders and Grevesse (1989). The dependence 
of the Fe abundance on the temperature of the ICM is 
weak (Fukazawa et al. 1998). ASCA revealed large-scale 
abundance gradients from AWM 7 and the Perseus clus- 
ter (Ezawa et al. 1997; 2001). De Grandi et al. (2004) de- 
rived the Fe abundances of nearby hot clusters observed 
with Beppo-SAX and found that the abundance profiles are 
strongly peaked for cool core clusters, whereas they remain 
constant for other systems. XMM-Newton and Chandra ob- 
servations show a spatial distribution of the Fe abundance 
of up to 0.3 ~ 0.4ri8o (e.g., Tamura et al. 2004; Vikhlinin et 
al. 2005; Baldi et al. 2007; Maughan et al. 2008; Leccardi 
& Molcndi 2008). At the center of most relaxed clusters, 
the Fe abundance decreases steeply outward. Outside the 
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central region, the Fe abundance decreases more gradually 
toward the outer regions. 

The Suzaku satellite first measured the Fe abundance 
of the ICM beyond O.SrigQ. (Fujita et al. 2008; Tawa 2008). 
Within 0.3-0.5rigo, Suzaku can derive the Fe abundances 
more accurately than XMM. The Fe abundance gradually 
decreases from the center to ^0.7rigo- 

In this paper, we describe our study of the Fe abun- 
dance in the ICM up to 0.3 ^ O.Srigo of 28 brightest clus- 
ters of galaxies observed with XMM-Newton. Some clusters 
have cool cores at their center (e.g. Makishima et al. 2001; 
Bbhringer et al. 2002). We must be careful when deriving el- 
emental abundances with a multi-temperature plasma, be- 
cause emission lines and a continuum spectrum are different 
functions of temperature. Therefore, we directly derived the 
strength of the Ka lines of Fe, considering the temperature 
dependence of the flux ratios of the Fe lines and the con- 
tinuum, and derived Fe abundances. We adopted an energy 
range in which the ICM dominates the background. 

In Section 2 we summarize the observations and data 
preparation. Section 3 describes our analysis of the data, 
and in Section 4 the Fe abundances are determined. We 
discuss our results in Secion 5. 

We adopt the new solar abundances given by Fodders 
(2003), according to which the solar Fe abundance with re- 
gard to H is 2.95xl0~^by number. This value differs from 
the photospheric value (4.68 xlO~^) given by Anders and 
Grevesse (1989). Considering the difference in He abun- 
dance between the two solar abundance tables, the Fe abun- 
dance yielded by the former is 1.5 times higher than that 
of the latter. We use Hq = 70 km/s/Mpc. Unless otherwise 
specified, errors are quoted with 68% confidence. 
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2. Target selection and observations 

In the XMM-Newton archival data, we selected the 28 
brightest clusters of galaxies with z<0.08 and total expo- 
sures above 7 ks after filtering out background flares. The 
samples are listed in Table 1. The ICM temperatures in the 
"the X-ray Brightest Abell Clusters (XBACs)" by Ebeling 
et al. (1996) range from 2.4 keV to 8.7 keV. Seventeen clus- 
ters show relatively relaxed morphologies with a cD galaxy 
at the X-ray peak. Hereafter, these are called cD clusters. 
Eleven clusters are merging clusters or relaxed clusters with 
two giant galaxies in the central region such as the Coma 
cluster. Hereafter, these clusters are called non-cD clusters. 

The virial radius of each cluster, rigo — 
1.95/1^-00 Vfc <T> /lOkeV, was calculated following 
Markevitch et al. (1998), by adopting an average tempera- 
ture < T > of that at 0.06-0.3ri8o from the X-ray peak of 
each cluster. 

We used the PN, MOSl, and M0S2 detectors to de- 
rive the Fe abundance of the ICM. We selected events with 
patterns smaller than 5 and 13 for the PN and MOS, re- 
spectively. 

Spectra were accumulated in concentric annular regions 
of 0-0.03ri8o, 0.03-0.06ri8o, 0.06-0.1ri8o, 0.1-0. 2ri8o, 0.2- 
0.3ri8o, and 0.3-0.5ri8o centered on the X-ray peak of 
each cluster. Here, the X-ray peaks were derived using the 
ewavelet tool of SAS-v8.0.0, and luminous point sources 
were excluded. The calculated X-ray peak of each cluster 
and the accumulation regions for the spectra are summa- 
rized as X-ray images in Figure [T] The spectra from MOSl 
and M0S2 were added. The background spectrum was cal- 
culated by integrating blank-sky data in the same detector 
regions. From deep-sky observations with the XMM, we se- 
lected the data with the background most similar to that 
of each cluster, and the faintest Galactic emission, after 
screening out background flare events in the data and the 
background, following Katayama et al. (2004). Then, the 
background was scaled using the count rate between 10 and 
12 keV. Most scaling factors are smaller than 10%, because 
we selected a blank-sky observation to match each cluster 
observation. 

The response matrix file and the auxiliary response file 
(ARF) corresponding to each spectrum were calculated us- 
ing SAS-vS.O.O. We used the XSPEC_vll.3.2ag package for 
our spectral analysis. 

3. Analysis 

3.1. Fe line strength 

In this study, the Fe abundances of the ICM are derived 
from the ratio of the flux in units of photons •cm~^s~^ of the 
Ka lines of He-like Fe (hereafter FfoHok) or H-like Fe (here- 
after FfcHk) to that of the continuum at 3.5-6.0 keV (here- 
after i^3.5-6)i since the systematic uncertainty in the con- 
tinuum flux in this energy band owing to the background is 
smaller than that around the Fe-K lines. Another advantage 
is that the dependence of the ratios on the plasma temper- 
ature is relatively weak within a certain temperature range. 
The temperature dependence of these ratios can constrain 
on the Fe abundances of multi-temperature plasmas. 

To obtain suitable statistics, we simultaneously fitted 
the raw annular spectra of a cluster and a deep field within 
an energy range of 5.0-7.2 keV with the sum of the ICM and 



background emission. The ICM component consists of ther- 
mal bremsstrahlung and two Gaussians for the Ka lines of 
He-like and H-like Fe. We modeled the background emission 
with a power-law model for the cosmic X-ray background 
(CXB), a "powerlaw/b" model for non X-ray background 
(NXB), and three Gaussians at 5.4 keV, 5.9 keV, and 6.4 
keV, which are the Ka lines of neutral Cr, Mn, and Fe, re- 
spectively. The "powerlaw/b" model is not folded through 
the ARF, and differs from a power-law. We obtained accept- 
able fits with a reduced ~ 1- Figure [5] shows representa- 
tive spectra of Abell 426 (A426, the Perseus cluster) fitted 
in this way. These spectra are fitted well, and the Ka line of 
He-like Fe is clearly detected even at 0.3-0.5ri8o. The back- 
ground line at 6.4 keV may cause a systematic uncertainty 
in the strength of the He-like Fe line in the MOS spectra at 
low-brightness regions. In contrast, the background lines in 
the PN spectra in this energy range are much weaker than 
those in the MOS. 

Despite the 6.4 keV background line, the PN and 
MOS gave mostly consistent ratios of FfoHck/^3.5-6 and 
^FeHK/^3.5-6- Hereafter we used the weighted average of 
these ratios from the PN and MOS to derive Fe abundances 
with relative errors smaller than 30%. 

3.2. Spectral fit with single-temperature model 

To convert the observed line flux to the Fe abundance, we 
derived the ICM temperatures. The first step in this anal- 
ysis was to fit the background-subtracted annular (pro- 
jected) spectra from the PN and MOS with a single- 
temperature vAPEC model (ISmith et al. 2001]) with pho- 
toelectric absorption fixed at the Galactic value (hereafter 
the IT model). The energy ranges of 1.4-1.6 keV for the 
MOS and 7.2-9.0 keV for the PN were discarded owing 
to strong instrumental lines. To account for the remaining 
background, we added a "powerlaw/b" model (Zhang et 
al. 2009) and fitted the spectra within the energy range of 
1.2-10.0 keV. Here, we fixed the photon index of the "pow- 
erlaw/b" component at 0.15 and 0.35 for the MOS and PN, 
respectively (Zhang et al. 2009). The abundances of C, N, 
and Al were fixed at 1 solar and those of O and Ne were 
assumed to have the same values as Mg. The abundances 
of Mg, Si, S, Ar, Ca, Fe, and Ni were allowed to vary. We 
did not use a more detailed modeling of the background 
(Snowden et al. 2008; Leccardi & Molendi 2008), since we 
mainly used an energy range over which the ICM emission 
is dominant. The ICM temperatures derived from the PN 
and MOS are consistent with each other within the order 
of several percent. We adopted the weighted average of the 
temperatures derived from the PN and MOS, as summa- 
rized in Tabled 

To derive the Fe abundance in each annular region, we 
included a 10% systematic uncertainty in temperatures, 
considering the uncertainties described below. When the 
spectra are fitted within the energy range of 0.8-10.0 keV, 
which includes Fe-L emission, the derived temperatures de- 
creased systematically by ^10%. We also fitted the spectra 
without the "powerlaw/b" model. Then, outside 0.2ri8o, 
the derived ICM temperature increased by 10~20% and 
the discrepancy between the PN and MOS became in- 
creased. Similarly, we also analyzed the same XMM data 
in the same way with SAS v7.0.0. In this case, when the 
derived temperatures are higher than several keV, the MOS 
systematically gives ~10% higher than those given by the 
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Fig. 1. Combined MOS and PN images of the target clusters. Exposures are not corrected and therefore the exposures 
of each field can be visualized. Circles have radii of 0.03, 0.06, 0.1, 0.2, 0.3, and O.Srigo- 
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Fig. 2. Representative spectra of AbcU 426 (the Perseus 
cluster) (black) and the adopted background (gray) at 0.3- 
0.5ri8o from MOS and PN fitted with Gaussians and con- 
tinuum. Best-fit Gaussians corresponding to the Ka-lines 
of He-like and H-like Fe are plotted as bold lines: those for 
background lines are plotted as thin dashed lines. 



PN. Systematic differences of ^10% in cluster temperature 
among Chandra, the PN, and the MOS were reported in 
relatively high-temperature clusters (Snowden et al. 2008). 
Reese et al. (2010) found that the ICM temperature de- 
rived from Chandra data changes systematically by ^^10 % 
between calibrations. 

3.3. Spectral fit with multi-temperature model 

To study the effect of the temperature structure of the ICM, 
we also fitted the background-subtracted spectra of the 
MOS at 0.5-10.0 keV, using a multi-temperature vAPEC 
model with photoelectric absorption (multi-T model), the 
Galactic emission, and the "powerlaw/b" model with a 
fixed photon index, as described in Section 3.2. The 
Galactic emission, which includes the local hot bubble, the 
Milky Way halo, and solar wind charge exchange, is em- 



pirically fitted with a two-temperature APEC model with 
redshift = (e.g. Lumb et al. 2002; Yoshino et al. 2009 ). 
Therefore, we used the two-temperature APEC model for 
the Galactic emission. The temperature of the one com- 
ponent was fixed at 0.1 keV, and that of the other was 
restricted within 0.2-0.4 keV. Here, we used only the MOS 
detector, because it has a better energy resolution and is 
more sensitive to temperature structure than the PN detec- 
tor. The multi-T model is a sum of 13 temperature compo- 
nents, because the spectrum of an isothermal plasma can 
be reproduced well by the sum of the spectra of two neigh- 
boring temperature plasmas with similar elemental abun- 
dances. The ratio of the temperatures of the two neighbor- 
ing components is fixed at 0.8. The hottest temperature 
of the multi-T model is fixed at the highest temperature 
derived using the IT model derived from the MOS within 
0.2ri8o. The elemental abundances of each element for the 
multi-T components are assumed to have the same values. 
We added a 3% systematic error to the count rate of each 
channel in the spectra. Then, the reduced reaches the 
reasonable values of ~1. 



4. Results 

4.1. Temperature profiles of the ICM 

The radial temperature profiles of the ICM are summarized 
in Figure|31 From 0.06ri8o to 0.2~0.3 riso, the temperature 
gradient of the ICM is fiat, whereas the temperatures of 
several cD clusters are lower within O.OGrigo than those at 
0.06~0.2ri8o, owing to the presence of cool cores. In con- 
trast, most non-cD clusters have flatter radial temperature 
profiles within 0.2~ O.Snso- Beyond 0.2~0.3ri8o, the ICM 
temperatures tend to decrease outward. These temperature 
profiles are consistent with previous Chandra and XMM- 
Newton measurements (e.g., Snowden et al. 2008; Vikhlinin 
et al. 2005). 

4.2. Fe abundance of the ICM 

Figure m shows the observed values of FfcHck/-F3.5_6 and 
^FcHk/-F3.5-6: plotted against the ICM temperature de- 
rived from the IT model. Within 0.03ri8o, ^FcHck/-F3.5-6 
is scattered from 0.03 to 0.2, with mild temperature 
dependence. Beyond 0.03ri8o, both FFcHeK/-F3.5-6 and 
-FFeHK/-F3.5-6 become closer to the theoretical relationships 
predicted by the APEC model assuming that the Fe abun- 
dance=0.5 solar. Beyond 0.3ri8o, the data for only several 
clusters have adequate statistics: the i^FcHcK/-p3.5_6 values 
of A426 (the Perseus cluster) and a non-cD cluster are close 
to the relationship for 0.5 solar, whereas those of two non- 
cD clusters are significantly lower. 

The temperature dependence of the theoretical ratio 
-FFeHeK/^3.5-6 for a given Fe abundance is fairly weak, 
within 20% for 2-6 keV, and starts to decrease above 6 
keV (Figure HJ. In contrast, the ratio of ^fcHk/^3.5-6 is 
nearly constant within 20% for 6.7-17 keV. This tempera- 
ture independence minimizes the systematic uncertainty in 
the Fe abundance because of uncertainties in the tempera- 
ture structure. 

To derive the Fe abundance, we converted the observed 
ratio of FfcHok/^3.5-6 using the theoretical ratios calcu- 
lated by the APEC model and the ICM temperature de- 
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Fig. 3. Radial profiles of the ICM temperature of cD clus- 
ters (upper panel) and non-cD clusters (lower panel). 
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The results are summarized in 



4.3. Radial profiles of Fe abundances 

The radial profiles of the Fe abundances in the ICM are 
shown in Figure [SJ Within O.OSrigo, the derived Fe abun- 
dances of the cD clusters are 0.7-1.6 solar. Non-cD clusters 
with large deviations from spherical symmetry, A1367 and 
A3562, also have high Fe abundances of ^ 1.5 solar. In con- 
trast, the other non-cD clusters have lower central values of 
Fe. Outside O.lrigo, both the cD and non-cD clusters have 
similar radial profiles of the Fe abundance. At 0.1-0.2rigo, 
the derived Fe abundances are '^0.5 solar. For the region at 
O.S-O.Sriso, although sufficient statistics are available for 
only six clusters, the derived Fe abundances are around 
0.3-0.5 solar. 

Figure [S] shows the ratios of the Fe abundances of 
each cluster at each annular region to that at 0.06-0. IrigQ. 
Except for the Centaurus cluster (A3526), the Fe abun- 
dances of the cD clusters have similar slopes: the Fe abun- 
dances within 0.03rigo are higher than those at 0.06-0. Irigo 
by a factor of 1.5. Beyond 0.06rigo, the radial profiles be- 
come flatter. 

The Fe abundances of three morphologically distorted 
clusters, A1367, A3558, and A3562, within O.OSrigo are 



Fig. 4. FfcHck/^3.5-6 (upper panel; circles) and 
-FFeHK/-F3.5-6 (lower panel; triangles) at 0-0.03rigo (red), 
0.03-0.06rigo (orange), 0.06-0. Irigo (magenta) 0.1-0.2rigo 
(green), 0.2-0.3rigo (blue), and 0.3-0.5rigo (black), plotted 
against the ICM temperature. Closed and open symbols in- 
dicate cD clusters and non-cD clusters, respectively. Black 
solid and dotted lines correspond to the theoretical ratios 
from the APEC model of JFeHeK/-F3.5-6 and FpeRK/ Fs.^-e, 
respectively, assuming Fe abundance = 0.3, 0.5, 1.0, 1.5, 
and 2.0 solar. Here, abundance ratios of elements are as- 
sumed to have the solar ratios by Lodders (2003). Those of 
the MEKAL model assuming Fe abundance = 0.5 solar are 
plotted as gray lines. 



higher than those at 0.06-0. Irigo by a factor of 1.5-2. The 
other non-cD clusters have similar Fe abundances within 
0.03rigo and 0.06-0. Irigo. 

4.4. Temperature dependence of the Fe abundance 

In Figure [7] the derived Fe abundances are plotted against 
the ICM temperature. There is no significant dependence 
on the plasma temperature, or between the cD and non-cD 
clusters. The scatter in the Fe abundance tends to be larger 
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Fig. 6. Ratios of Fe abundances derived from 
^FcHok/^3.5-6 at 0-0.03ri8o (red), 0.03-0.06ri8o (orange), 
0.1-0. 2ri8o (green), 0.2-0.3ri8o (blue), and 0.3-0.5ri8o 
(black) to those at 0.06-0. lri8o, plotted against the ICM 
temperature at 0.06-0. lri8o. Upper and lower panels 
correspond to the cD and non-cD clusters, respectively. 



Fig. 5. Radial profiles of the Fe abundances of cD clusters 
(upper panel) and non-cD clusters (lower panel). 

Table 3. Averages of the relative differences in the de- 
rived Fe abundances [(Fe2 — Fei)/Fei] and their root- mean- 
square scatters. Cool (hot) samples have local ICM temper- 
atures below (above) 5 keV. 



Fei 


FC2 


sample 


average 


scatter 
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0.10 


FeHK, IT 


FeHK, multi-T 


hot 


-0.04 


0.07 



in lower temperature clusters. When the ICM temperature 
is higher than 5 keV, the scatter in the absolute values of 
the Fe abundances at 0.06-0.3ri8o is small. 

4.5. Systematic uncertainties in Fe abundance 

In determining the abundance, we must examine uncertain- 
ties such as the dependence on the temperature structure 
and plasma codes. 

When the ICM temperature derived from the IT model 
was higher than 5 keV, the observed -FfcHk/^3.5-6 was con- 
verted to the Fe abundance as described in Section 4.2. 
The derived Fe abundances agree well, within statistical er- 
rors, with those derived from FFoHeK/^3.5-6 (upper panel 
in Figure |S1 Table 2). On average, the difference between 
the Fe abundances derived from the two lines is only a few 
%. At the center of the cool core, where the temperature 
and abundance gradients of the ICM are steep, the Fe abun- 
dances derived from the He-like line can be expected to 



be higher than those from the H-like line. However, outside 
the cool core, where the temperature and abundance gra- 
dients weaken, it is reasonable that He-like and H-like lines 
yield the same Fe abundances. Because the temperature de- 
pendence of the He-like and H-like Fe lines are completely 
different, the consistency of the Fe abundances outside the 
cool cores indicates a small systematic uncertainty owing 
to the uncertainties in the ICM temperature. 

To examine the systematic effect of the assumptions 
about the temperature structure, i.e., the IT model or the 
multi-T model, we derived the Fe abundance from the ob- 
served -FFeHeK/-F3.5-6 and FfoHk/-F3.5-6 assuming the best- 
fit multi-T model. For both He-like and H-like lines, the 
multi-T model yields mostly the same Fe abundances as 
the IT model, as summarized in Figure [9] and Table [21 
In Figure [HI the error bars arise from a systematic uncer- 
tainty of 10% in the ICM temperature, because in most 
cases these errors are much larger than those from the un- 
certainty in the multi-T model. Statistical errors in the Fe 
line strengths are not considered here, because we use the 
same line strength. As summarized in Table [3l the differ- 
ences in the Fe abundances from -FpeHeK/^a.s-e are smaller 
than several percent particularly below 5 keV owing to the 
weaker temperature dependence. Above 5 keV, the system- 
atic differences are still smaller than ~10%. The differences 
in the Fe abundances from FFeHK/-?^3.5-6 are also about sev- 
eral percent owing to the weaker temperature dependence. 

Figure [5] also compares the Fe abundances from the two 
lines, assuming the best-fit multi-temperature model. Here, 
statistical errors in the line strengths and a 10% systematic 
uncertainty in the temperature are considered. In most clus- 
ters, the two lines yield the same Fe abundances, reflecting 
that the IT and multi-T models yielded the same Fe abun- 
dances. 
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Fig. 7. Fe abundances derived from i^FcHcK/^3.5-6 using 
the IT APEC model, plotted against the loeal ICM temper- 
ature. The indieations of symbols and colors are the same 
as in Figure m 



Higher-temperature components may exist in the ICM. 
Therefore, we left the highest temperature of the multi- 
T model unconstrained and fitted the spectra from the 
PN and MOS. For the ICM temperatures from the IT 
model below 5 keV we obtained almost the same Fe abun- 
dances within several percent. When the ICM temperature 
is above 5 keV, the Fe abundances using the best-fit multi-T 
model from the PN spectra were unchanged. However, us- 
ing the MOS multi-T model, the Fe abundances from the 
He-like line increased by a several tens of percent, whereas 
those from the H-like line remained the same within 10%. 
This discrepancy arises from a systematic uncertainty of 
several percent in the response matrices of the PN and 
MOS. Considering the weak temperature dependence of the 
^FeHK/-p3.5-6 at higher temperatures, the Fe abundances 
from the H-like line are more reliable outside the cool core 
regions. 

Rasia et al. (2008) studied bias or systematic effects in 
the X-ray spectra using the spectra simulated by the X-ray 
Map Simulator and found that the Fe abundance is recov- 
ered with high accuracy for hot (T > 3 keV) and cold (T < 

2 keV) systems. At intermediate temperatures (2 keV< T < 

3 keV), they found systematic overestimates. The problem 
occurs because of a transition between the relative impor- 
tance of the Fe-L and Fe-K lines. Because we here derived 
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Fig. 8. Fe abundances derived from FpeHK/^a.s-e ver- 
sus those from FfoHck/^3.5-6j using the IT model (up- 
per panel) and the best-fit multi-T model (lower panel). A 
systematic uncertainty of 10% in the temperature is con- 
sidered. The indications of colors are the same as in Figure 
4. Closed and open circles correspond to cD and non-cD 
clusters, respectively. 



the Fe abundances from the Fe-K lines, the observed values 
do not include those of the coolest components (below ~ 1.7 
keV). The fraction of the emission measure for temperature 
components below this temperature is smaller than several 
percent except in the innermost region of the coolest cool 
cores of the Virgo cluster and A 262, where the fractions 
are 30%. 

To examine the systematic effect of the plasma code on 
abundance determination, we plotted the temperature de- 
pendence of the ratio of line strength to i^a.s-e in Figured 
using the MEKAL model (Mewe et al. 1995, 1996; Kaastra 
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kT (keV) 



Fig. 9. Differences between the derived Fe abundances us- 
ing the best-fit multi-T model and those using the IT 
model, plotted against the ICM temperature. Error bars 
are the result of uncertainties in the ICM temperature de- 
rived from the IT model, including systematic uncertainties 
of 10%. Upper and lower panels correspond to Fe abun- 
dances derived from FfcHck/-F3.5-6 and FfcHk/-F3.5-6, re- 
spectively. The indications of the colors are the same as in 
Figure 4. 
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Fig. 10. Fe abundance profiles of A426 (green), A1656 
(blue), and the averages of those of cD clusters (red 
solid diamonds) and non-cD clusters (magenta dashed di- 
amonds) in our sample, compared with the average profile 
of z=0.1^0.3 clusters (black closed circles) from Leccardi 
& Molendi (2008). Black dashed line corresponds to the Fe 
abundance derived from the numerical simulations without 
AGN feedback by Fabjan et al. (2008). 



clusters have systematically lower Fe abundances within 
O.lrigo, and the values similar to those of cD clusters be- 
yond this radius. The Coma cluster also has a fiatter radial 
profile within 0.2ri8o, but at a lower value of 0.4 so- 
lar, and beyond this radius, the Fe abundance gradually 
decreases to ~ 0.3 solar. 



1992; Liedahl et al. 1995). The difference between the 
MEKAL and APEC models is less than 10%. 

Non-solar abundance ratios of He, C, N, and O to Fe 
cause a bias in the derived Fe abundance, since these ele- 
ments change the continuum level. The difference in the He 
abundance between the solar abundance tables by Lodders 
(2003) and Anders & Grevesse (1989) yields a bias of sev- 
eral percent. The bias resulting from a possible high He- 
abundance in the ICM is discussed in detail in Ettori & 
Fabian (2006) and Bohringer & Werner (2010). 

In summary, in most cases, the systematic uncertainty 
in the derived Fe abundances caused by the uncertainties in 
the temperature structures and atomic data is about 10%. 

4.6. Weighted average of Fe abundance profiles 

We calculated the weighted average of the Fe abundances 
derived from FFcHcK/-f^3.5-6 and i^FcHK/^3.5-6, using the 
IT APEC model within the same radial region in the 
units of riso . We divided the clusters except for A426 (the 
Perseus cluster) and A1656 (the Coma cluster) according 
to whether they are cD or non-cD clusters. 

The results are summarized in Table |4] and Figure 1101 
The radial profile of the average Fe abundance of the cD 
clusters becomes flatter from 0.1 to 0.3 rigp at 0.4-0.5 so- 
lar. The Perseus cluster has a 10% smaller radial profile 
than the other cD clusters within 0.3ri8o, and a fiatter Fe 
abundance profile up to 0.5ri8o. In contrast, the non-cD 



Table 4. Weighted average of Fe abundances derived from 
-FfoHok/-F3.5_6 and i^FoHK/^3.5-6, in the units of the solar 
abundance. 



radius" 



cDs" 



A 426 



non-CDs'^ 



A 1656 



0.00-0.03 
0.03-0.06 
0.06-0.10 
0.10-0.20 
0.20-0.30 



0.88±0.02 
0.73±0.01 
0.60±0.01 
0.48±0.01 
0.42±0.04 



0.78±0.01 
0.66±0.05 
0.55±0.05 
0.44±0.04 
0.41±0.05 



0.30-0.50 0.35±0.10 0.48±0.06 



0.69±0.02 
0.57±0.02 
0.52±0.02 
0.50±0.02 
0.44±0.03 
0.32±0.04 



0.45±0.03 
0.41±0.03 
0.39±0.03 
0.39±0.03 
0.31±0.03 
0.25±0.03 



in unit of rigo. 

the cD clusters except A 426 

the non-cD clusters except A 1656 



4. 7. Comparison witfi previous papers and redsfiift evolution 

In Figure [lU] we compare our results with those obtained by 
Leccardi & Molendi (2008). Here, the Fe abundances were 
scaled using Lodders (2003), considering the difference in 
the solar values of He and Fe. Leccardi and Molendi (2008) 
analyzed ^ 50 hot [kT > 3.3 keV) intermediate-redshift 
(z — 0.1 ^ 0.3) clusters of galaxies observed with XMM- 
Newton. Although we use samples different from those em- 
ployed by Leccardi & Molendi (2008), our mean profiles of 
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the cD clusters and those of Leccardi & Molendi (2008) have 
remarkably similar slopes, but the absolute values differ by 
10~20%. 

Maughan et al. (2008) analyzed redshift evolution up 
to z ~ 1 of the Fe abundance in the ICM using Chandra 
observations. The mean value of Maughan et al. (2008) at 
0.15-lr5oo of z < 0.5 clusters is ~ 0.5 solar using the solar 
abundance by Lodders (2003), and is consistent with the 
average of the cD clusters. 

These agreements indicate that the evolution of the Fe 
abundance with redshift from z = to z = 0.5 is consis- 
tent with no evolution within our systematic uncertainties. 
The expected evolution according to Ettori (2005), which 
is about 30% from 0.5 to z^^O, tends to be larger than 
our observations. 



5. Discussion 

5.1. Flatter radial Fe abundance profile beyond O.lriso 

Beyond O.lrigo, well outside the cool cores, both the cD and 
non-cD clusters have similar Fe abundance profiles with a 
relatively small scatter. Within O.l-O.Srigo, the Fe abun- 
dances of the sample clusters are 0.4~0.5 solar. There is no 
dependence on the ICM temperature, i.e., system richness. 
In other words, outside the central regions of clusters, the 
Fe abundance in the ICM may be universal in clusters of 
galaxies, and clusters of galaxies may have universal metal- 
enrichment histories. Beyond 0.3ri8o, the Fe abundances 
become 0.2-0.5 solar. These profiles smoothly connect with 
the Suzaku observations of the Fe distributions in the ICM 
beyond 0.5ri8o (Fujita et al. 2008; Tawa 2008). 

The distribution of Fe abundance in the ICM has been 
modeled via hydrodynamic simulations (Fabjan et al. 2008, 
Kapferer et al. 2007). Figure [TOl also compares our mean 
Fe abundance profiles with that derived by Fabjan et al. 
(2008). Beyond O.lrigo, the profile simulated by Fabjan et 
al. (2008) continues to decrease outward, whereas the ob- 
served profile becomes flatter. 

The iron-mass-to- light ratio (IMLR), i.e., the ratio of Fe 
mass in the ICM to total stellar luminosity, is important for 
studying the origin of Fe, because Fe is synthesized in stars. 
The integrated values of the IMLR profiles of the kT ^ 3 
keV clusters, Abell 262, and AWM7, observed with Suzaku 
(Sato et al. 2008, 2009), increase outward from O.It'iso to 
0.3-0.4ri8o. This result indicates that Fe in the ICM ex- 
tends farther than stars in the outer regions at least in 
these two clusters. 

A simple explanation is that a significant fraction of 
Fe is synthesized in an early phase of cluster evolution. If 
metal enrichment occurs after the formation of the clusters, 
the metal distribution is expected to follow the stellar dis- 
tribution, whereas if it occurs before cluster formation, we 
expect a flatter distribution. For example, simulations by 
Kapferer et al. (2007) yielded a flatter Fe profile from SN II. 
That the measurements of the Fe abundance excluding the 
central region by Maughan et al. (2008) are consistent with 
no evolution at least up to z=0.7 also supports the early 
synthesis of Fe. Feedback by active galactic nuclei (AGN) 
may also change the abundance profile. With AGN feed- 
back, a flatter radial profile beyond O.lrigo is also obtained 
in the simulation by Fabjan et al. (2010). 

A426 (the Perseus cluster) shows a completely flat Fe 
profile from O.lrigo to 0.5rigo: at 0.3-0.5rigo, the Fe abun- 



dance is higher than in the other clusters. Within the an- 
nular region of 0.3-0.5ri8o, the observed fields of view do 
not cover the entire annular region (Figure 1) and lack the 
north and south regions. Several clusters, including A426, 
form the Perseus supercluster, which is elongated in an east- 
west direction (Gregory et al. 1981). The observed very flat 
Fe abundance may reflect the history of cluster formation 
from filaments of the large-scale structure of the Universe. 

5.2. Fe profiles within the cool core regions 

Fe in the gas in the cool core of clusters is a mixture of 
that present in the ICM and later supplies from their cD 
galaxies. The latter contains Fe synthesized by SN la and 
originating from stars through stellar mass loss, because Fe 
is synthesized by both SN la and SN II. According to abun- 
dance patterns from XMM-Newton observations of the cool 
cores (Matsushita et al. 2003, 2007; Werner et al. 2006; de 
Plaa et al. 2006; Simionescu et al. 2008, 2009) and Suzaku 
observations of the cD galaxies of AWM7 and A262 (Sato 
et al. 2008, 2009), the number ratio of SNe II to SNe la 
is estimated to be 3~4, and the Fe in the cluster core was 
synthesized mostly by SN la. 

Within cool core regions, clusters have negative abun- 
dance gradients (e.g. Fukazawa et al. 2000, de Grandi et al. 
2004). Rossetti & Molendi (2010) and Leccardi, Rossetti, 
& Molendi (2010) also found central metal abundance ex- 
cess in low-entropy core systems. Within these regions, 
observations indicate that metals ejected from cD galaxies 
have a more extensive distribution compared with stars. 
Analyzing eight cool-core clusters observed with XMM, 
David & Nulsen (2008) found that Fe in the ICM is more 
extended than the stellar mass, and excess Fe within the 
central 100 kpc could have been produced by type la super- 
nova in cD galaxies over the past 3-7 Gyr. The IMLR val- 
ues in AWM7 and A262 observed with Suzaku have much 
steeper positive gradients within the cool cores than those 
beyond 0.1ri8o (Sato et al. 2008; 2009). 

Processes such as jets from a central AGN as described 
below, or the "sloshing" of cD galaxies in the cluster's grav- 
itational potential (e.g. Markevitch et al. 2001), may eject 
metals from cD galaxies. On the basis of 2D temperature 
and metallicity maps of M 87 (Simionescu et al. 2008) and 
Hydra A (Simionescu et al. 2009), uplift of Fe by the AGN 
was discussed in detail. David & Nulsen (2008) found that 
both turbulent diffusion of entropy and dissipation are im- 
portant heating mechanisms in cluster cores. Rebusco et 
al. (2006) compared the observed Fe abundance profiles 
with the predictions of a model involving metal ejection 
from the brightest galaxy and the subsequent diffusion of 
metals by stochastic gas motions. Roediger et al. (2007) 
studied the effect of bubble-induced motions on metallicity 
profiles through 3D hydro simulations. The new simulation 
by Fabjan et al. (2010) showed a less steep central Fe peak 
including feedback from AGN in cD galaxies. The observed 
small scatter in the Fe abundance slope within 0.1ri8o in the 
cD clusters suggests a common process of metal supply and 
ejection in the cD galaxies. The next Japanese X-ray satel- 
lite, ASTRO-H, will directly reveal line broadening caused 
by turbulence in the cores of clusters. 

Peng & Nagai (2009) pointed out that the He/H abun- 
dance ratio could be four times the solar abundance at the 
center of massive clusters. Such a high He abundance causes 
an underestimation of the Fe abundance by a factor of two 
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(Ettori & Fabian 2006). However, in actual clusters, metals 
at the center of the cluster come from mass loss by central 
galaxies, and in addition, more extensively distributed com- 
pared with stars. Therefore, the centers of clusters should 
not be static, as assumed in the calculations of Peng & 
Nagai (2009). 

5.3. Mixing of the ICM in non-cD clusters 

Most of clusters without cool cores tend to show no central 
abundance excess as found in Fukazawa et al. (2000); de 
Grandi et al. (2004). Leccardi, Rossatti & Molendi (2010) 
found that some medium- and high-entropy core systems 
have central abundance excess, most of which are undergo- 
ing a phase of rapid dynamical changes, and suggested that 
these clusters originate from low-entropy core systems. 

Our results are consistent with the discussion in 
Leccardi, Rossatti & Molendi (2010). The non-cD clusters 
in this paper include merging clusters in various stages, 
from heavily distorted morphologies to relatively relaxed 
ones without cool cores. During cluster merging, mixing of 
the ICM could destroy the central Fe peak. Clusters with 
the highest central Fe abundances have highly disturbed 
morphologies. In the first stage of merging, the Fe peak 
and the cool core may remain. Then, mixing of the ICM 
destroys the Fe peak, and the Fe abundance becomes flat 
at the same level as in cD clusters at 0.1-0.5rigo. 

6. Summary and conclusion 

We derived radial profiles of the Fe abundance of the ICM 
in nearby {z < 0.08) clusters observed with XMM-Newton. 
The Fe abundances of the ICM were derived from the ra- 
tio of the flux of the Ka lines of He-like or H-like Fe to 
those of the continuum at 3.5-6.0 keV, because the system- 
atic uncertainty in the continuum flux in this energy band 
owing to the background is smaller. The temperature de- 
pendence of these ratios constrains the Fe abundances of 
multi-temperature plasmas. 

In cluster core regions (< O.lrigo), the observed Fe 
abundances of cD clusters show similar radial profiles. The 
less-peaked abundance profile compared with the light in 
the central region indicates the ejection of metals from cD 
galaxies. 

In the outer regions, 0.1-0.2ri8o, Fe abundances of 0.4- 
0.5 solar appear to be universal with no temperature de- 
pendence. The observed flatter Fe abundance proflles of the 
cD clusters beyond O.lrigg indicate early metal enrichment. 

Chandra and XMM cannot reach beyond 0.5rigo ow- 
ing to a high, unstable particle background. Information is 
available for only about 10% of the cluster volume, and the 
majority has not yet been revealed. The total amount of 
Fe synthesized in galaxies can be derived only by the pre- 
cise abundance measurements beyond 0.5rigo. Because of 
its low background, Suzaku is the only satellite available 
for the next several years to study clusters of galaxies up 
to the virial radius. With Suzaku, Fe abundances will be 
determined out to a radius of 0.7^ O.SrigQ. 
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Table 1. Cluster samples in the XMM- Newton archival data 





z < T >" 


N'' 


type 


obsid (exposures of PN, MOSl and M0S2)'' 






(keV) 


(10 


cm J 


(ks, ks, ks) 




A262 


0.016 


2.5 


5.4 


cU 


0109980101 


(19.0 23.6 23.9) 0504780201 (24.9 31.6 33.1) 


Virgo 


0.004 


2.8 


2.6 


CU 


0114120101 


(30.0 37.7 38.5) 0200920101 (69.6 76.9 78.5) 












0106060101 


(5.0 9.2 9.2) 0106060201 (1.6 5.5 5.0) 0106060301 (5.4 7.1 7.9) 












0141570101 


(17.4 23.2 24.7) 0112550701 (1.0 0.2) 0112552101 (7.6 13.4 13.5) 












0106060401 


(7.3 10.7 11.1) 0145800101 (5.9 37.8 0.2 0.2) 0106060501 (10.9 15.3 15.4) 












0200650101 


(42.7 53.2 53.5) 0112840101 (14.6 18.1 17.6) 0106060601 (8.6 10.5 10.7) 












0112550801 


(2.9 13.2 12.1) 0106060701 (1.2 12.4 5.9) 0108860101 (16.7 20.6 19.2) 












0106061401 


(4.9 8.0 8.3) 


A4Uo8 


0.029 


3.0 


i.o 


cU 


0204460101 


(27.8 29.3 28.9) 


AiUOU 


0.012 


3.0 




non-cD 


0206230101 


(28.5 37.9 38.9) 


A onKO 
A/Uo/ 


0.035 


3.1 


3.0 


CU 


0109920101 


(26.4 28.9 29.0) 


A1367 


0.021 


3.3 


2.4 


non-cD 


0061740101 


(23.7 31.8 31.4) 0005210101 (20.8 27.8 28.8) 


A ^on 


0.057 


3.5 


A ^ 

4. / 


CU 


0109980301 


(14.0 19.1 22.2) 


A2589 


0.042 


3.6 


3.9 


cU 


0204180101 


(21.8 26.3 27.5) 


AW M / 


0.017 


3.6 


y.o 


CU 


0135950301 


(27.5 30.7 30.9) 


MR W OS 


0.045 


3.7 


3.0 


CU 


0109930101 


(0.0 37.7 38.0) 


A on£?o 
AzUbo 


0.036 


3.9 


3.0 


CU 


0200120401 


(5.0 8.2 5.0) 


A3526 


0.010 


4.0 


8.2 


cU 


0046340101 


(35.2 44.9 43.5) 0406200101 (90.5 111.7 111.7) 0504360101 (26.4 28.2 30.9) 












0504360201 


(32.9 34.0 34.6) 


A2199 


0.030 


4.2 


0.9 


cU 


0008030601 


(0.5 6.9 4.7) 0008030301 (2.9 4.7 4.7) 0008030201 (13.3 14.5 14.6) 


A A Ciii 

A49b 


0.033 


4.4 


A G 

4.6 


cU 


0135120201 


(11.8 16.2 15.9) 0506260301 (44.5 57.0 59.5) 0506260401 (44.8 55.8 58.6) 


A 1644 


0.048 


4.6 


4.8 


non-cD 


0010420201 


(12.7 14.1 14.1) 




0.050 


4.8 


A n 
4.U 


non-cD 


0105261301 


(37.7 39.0 39.5) 0105261501 (4.8 22.4 22.4) 0105261801 (4.7 10.9 10.4) 












0105261701 


(0.0 20.6 20.4) 0105261601 (16.3 18.3 19.5) 0105261401 (9.0 12.8 12.4) 


A o K o 


0.048 


5.4 


A O 

4.2 


non-cD 


0107260101 


(38.2 43.1 42.5) 


A OCO^ 

A6o2l 


0.016 


5.5 


on o 


non-cD 


0204250101 


(2.4 6.2 5.8) 0208010101 (9.1 12.6 13.0) 0208010201 (10.9 13.6 14.3) 


Ai lyo 


0.062 


5.8 


1 o 

i.2 


cU 


0097820101 


(25.1 38.3 37.9) 0109070201 (53.0 54.5 54.6) 0205190101 (26.9 29.5 28.8) 












0205190201 


(20.8 22.6 23.2) 


A O K 
AOO 


0.052 


5.8 


3.1 


CU 


0065140101 


(9.8 12.4 12.2) 0065140201 (9.2 12.1 12.4) 


A 

AjddY 


0.053 


6.0 


4.6 


non-cD 


0105260101 


(0.8 7.3 6.8) 0206850101 (52.3 59.2 60.0) 0105260601 (17.7 23.1 23.3) 












0105260301 


(12.8 17.0 16.2) 0105260401 (11.5 16.3 16.5) 0105260501 (11.2 12.6 12.6) 












0105260201 


(15.1 19.0 18.4) 


A /I o^; 
A4ZD 


0.018 


6.1 


14.6 


cU 


0085110101 


(47.5 48.8 50.7) 0204720101 (11.9 15.1 15.2) 0085590201 (38.4 41.7 40.0) 


(Perseus) 










0305720101 


(10.0 13.0 12.6) 0305720301 (15.9 19.9 21.2) 0151560101 (18.0 24.9 25.1) 












0204720201 


(20.7 24.5 24.5) 0405410201 (9.4 17.3 25.5) 0405410101 (12.2 18.2 16.1) 


A2256 


0.058 


6.3 


4.1 


non-cD 


0141380101 


(6.4 9.2 8.7) 0141380201 (10.6 13.2 12.4 ) 












0112951501 


(5.9 10.6 10.9) 0112951601 (7.2 12.7 12.7) 0112950601 (8.4 12.3 11.8) 












0112500101 


(21.8 25.0 25.2) 




0.040 


6.5 


4.4 


CU 


0086950201 


(14.3 25.4 25.2) 


A onon 
A/U29 


0.077 


7.4 


3.2 


cU 


0111270201 


(9.3 12.1 12.4) 


AiDOD 


0.023 


7.8 


0.9 


non-cD 


0300530201 


(0.4 3.6 3.4 ) 0153750101 (19.2 21.3 21.8) 0124711401 (14.7 17.4 16.5) 


(Ooma) 










0300530101 


(18.6 20.9 21.0) 0300530301 (26.5 30.5 30.2) 0300530701 (16.4 24.4 24.6) 












0300530401 


(13.3 18.8 19.6) 0300530601 (16.7 22.6 22.1) 0300530501 (19.2 24.7 24.6) 












0124710901 


(16.7 21.6 22.1) 0124712001 (9.2 15.4 14.6) 0124710501 (10.7 24.4 24.8) 












0124710601 


(4.7 9.5 9.0) 0204040101 (69.7 73.0 74.0) 0204040301 (43.7 52.9 53.2) 












0204040201 


(54.8 61.8 64.8) 0124710801 (18.4 24.5 24.5) 0124710401 (8.2 5.9 5.6) 












0124712501 


(23.8 27.7 27.7) 0124712401 (11.0 17.4 18.1) 0124711101 (13.2 19.3 21.1) 












0124712201 


(18.2 26.7 26.2) 0124712101 (20.3 26.0 26.1) 0124710301 (10.6 13.4 14.0) 












0124710101 


(23.7 30.6 30.8) 


A3266 


0.055 


8.4 


1.5 


non-cD 


0105262201 


(3.4 0.0 3.6) 0105260901 (18.6 24.4 24.3) 0105260801 (16.2 19.6 20.4) 












0105262501 


(3.5 7.6 7.7) 0105262001 (2.9 7.6 7.6) 0105260701 (16.1 20.6 20.5) 












0105261101 


(8.0 12.7 11.1) 0105262101 (4.4 5.6 7.2) 


A754 


0.054 


8.6 


4.7 


non-cD 


0112950301 


(8.0 13.3 13.7) 0112950401 (7.6 14.3 14.3) 0136740201 (2.3 6.0 6.7) 












0136740101 


(12.5 14.2 14.2) 



°: Adopted average temperature of ICM used to calculate rigo- 

':The Galactic value of the hydrogen column density (Dickey & Lockman 1990) 

■^i "cD" corresponds to cD clusters, which are relaxed clusters with a cD galaxy at their center. Other clusters are classified as 
"non-cD." 

''lobsid (a unique XMM-Newton observation identifier) and exposure times of PN, MOSl, and MOS2, respectively, after screening 
out the background flares. 



K. Matsushita: Fe abundance in the Intracluster Medium 
Table 2. Results of spectral fitting of annular spectra. 
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cluster r (riso) fcT" (keV) Fe (solar) from FfcHck/^3.5-6 Fe (solar) from i^FcHK/-F3.5-6 
Best-fit" ± A^^t.t Af^, Best-fjt" ± A.^tat 



A262 


0, 


.00- 


-0, 


.03 


1, 


.95±0, 


.03 


1, 


.31±0, 


.28-0.26 




A262 


0, 


.03- 


-0, 


.06 


2, 


.36±0, 


.04 


0, 


.76±0, 


.09-0.03 




A262 


0, 


.06- 


-0, 


.10 


2, 


.43±0, 


.06 


0, 


.62±0, 


.08-0.02 




A262 


0, 


.10- 


-0, 


.20 


2, 


.12±0, 


.08 


0, 


.90±0, 


.21-0.13 




Virgo 


0, 


.00- 


-0, 


.03 


2, 


.11±0, 


.00 


0, 


.86±0, 


.09-0.11 




Virgo 


0, 


.03- 


-0, 


.06 


2, 


.47±0, 


.01 


0, 


.52±0, 


.05-0.02 




A4038 


0, 


.00- 


-0, 


.03 


3 


.12±0, 


.06 


0, 


.79±0, 


.06-hO.OO 




A4038 


0, 


.03- 


-0, 


.06 


3 


.22±0, 


.08 


0, 


.66±0, 


.05-^0.00 




A4038 


0, 


.06- 


-0, 


.10 


3 


.04±0, 


.07 


0, 


.57±0, 


.05+0.00 




A4038 


0, 


.10- 


-0, 


.20 


2, 


.81±0, 


.10 


0, 


.55±0, 


.08+0.01 




A1060 


0, 


.00- 


-0, 


.03 


3 


.35±0, 


.03 


0, 


.67±0, 


.03-0.01 




A1060 


0, 


.03- 


-0, 


.06 


3 


.27±0, 


.04 


0, 


.55±0, 


.02-0.00 




A1060 


0, 


.06- 


-0, 


.10 


2, 


.93±0, 


.05 


0, 


.52±0, 


.05-0.01 




A1060 


0, 


.10- 


-0, 


.15 


2, 


.79±0, 


.12 


0, 


.46±0, 


.10-0.01 




A2052 


0, 


.00- 


-0, 


.03 


2, 


.57±0, 


.04 


1, 


.28±0, 


.17-0.07 




A2052 


0, 


.03- 


-0, 


.06 


3 


.12±0, 


.06 


0, 


.84±0, 


.05-0.00 




A2052 


0, 


.06- 


-0, 


.10 


3 


.12±0, 


.07 


0, 


.81±0, 


.06-0.00 




A2052 


0, 


.10- 


-0, 


.20 


3 


.01±0, 


.07 


0, 


.69±0, 


.06-0.01 




A1367 


0, 


.00- 


-0, 


.03 


2, 


.83±0, 


.23 


1, 


.71±0, 


.45-0.15 




A1367 


0, 


.03- 


-0, 


.06 


3 


.37±0, 


.18 


0, 


.78±0, 


.11+0.00 




A1367 


0, 


.06- 


-0, 


.10 


3 


.33±0, 


.17 


0, 


.71±0, 


.09-0.01 




A1367 


0, 


.10- 


-0, 


.20 


3 


.36±0, 


.14 


0, 


.45±0, 


.07+0.00 




A1367 


0, 


.20- 


-0, 


.30 


3 


.22±0, 


.30 


0, 


.71±0, 


.17-0.02 




A780 


0, 


.00- 


-0, 


.03 


3 


.41±0, 


.09 


0, 


.72±0, 


.06+0.02 




A780 


0, 


.03- 


-0, 


.06 


3 


.67±0, 


.10 


0, 


.56±0, 


.05+0.01 




A780 


0, 


.06- 


-0, 


.10 


3 


.62±0, 


.12 


0, 


.41±0, 


.07+0.01 




A780 


0, 


.10- 


-0, 


.20 


3 


.93±0, 


.12 


0, 


.36±0, 


.04+0.00 




A780 


0, 


.20- 


-0, 


.30 


4, 


.33±0, 


.39 


0, 


.41±0, 


.13-0.02 




A2589 


0, 


.00- 


-0, 


.03 


3 


.48±0, 


.15 


1, 


.50±0, 


.17-0.01 




A2589 


0, 


.03- 


-0, 


.06 


3 


.59±0, 


.14 


1, 


.06±0, 


.10-0.01 




A2589 


0, 


.06- 


-0, 


.10 


3 


.50±0, 


.14 


0, 


.92±0, 


.10-0.01 




A2589 


0, 


.10- 


-0, 


.20 


3 


.34±0, 


.14 


0, 


.79±0, 


.10-0.01 




AWM7 


0, 


.00- 


-0, 


.03 


3 


.44±0, 


.05 


1, 


.12±0, 


.06+0.01 




AWM7 


0, 


.03- 


-0, 


.06 


3 


.78±0, 


.07 


0, 


.89±0, 


.04+0.00 




AWM7 


0, 


.06- 


-0, 


.10 


3 


.70±0, 


.08 


0, 


.62±0, 


.03+0.00 




AWM7 


0, 


.10- 


-0, 


.19 


3 


.67±0, 


.10 


0, 


.63±0, 


.04+0.00 




MKW3s 


0, 


.00- 


-0, 


.03 


3 


.18±0, 


.10 


1, 


.06±0, 


.11+0.01 




MKW3s 


0, 


.03- 


-0, 


.06 


3 


.67±0, 


.12 


0, 


.88±0, 


.08+0.00 




MKW3s 


0, 


.06- 


-0, 


.10 


3 


.78±0, 


.15 


0, 


.57±0, 


.06+0.00 




MKW3s 


0, 


.10- 


-0, 


.20 


3 


.47±0, 


.15 


0, 


.55±0, 


.08+0.01 




A2063 


0, 


.00- 


-0, 


.03 


3 


.29±0, 


.20 


1, 


.02±0, 


.32-0.01 




A2063 


0, 


.03- 


-0, 


.06 


3 


.66±0, 


.21 


1, 


.04±0, 


.20+0.01 




A2063 


0, 


.06- 


-0, 


.10 


4, 


.05±0, 


.39 


0, 


.61±0, 


.14-0.01 




A3526 


0, 


.00- 


-0, 


.03 


2, 


.58±0, 


.01 


1, 


.63±0, 


.18-0.14 




A3526 


0, 


.03- 


-0, 


.06 


3 


.90±0, 


.03 


0, 


.72±0, 


.02+0.00 




A3526 


0, 


.06- 


-0, 


.10 


3 


.81±0, 


.03 


0, 


.57±0, 


.03+0.01 




A3526 


0, 


.10- 


-0, 


.20 


3 


.19±0, 


.12 


0, 


.36±0, 


.03+0.01 




A2199 


0, 


.00- 


-0, 


.03 


3 


.63±0, 


.05 


0, 


.78±0, 


.04+0.04 




A2199 


0, 


.03- 


-0, 


.06 


4, 


.25±0, 


.08 


0, 


.60±0, 


.04-0.01 




A2199 


0, 


.06- 


-0, 


.10 


4, 


.28±0, 


.12 


0, 


.60±0, 


.04-0.01 




A2199 


0, 


.10- 


-0, 


.20 


4, 


.33±0, 


.11 


0, 


.45±0, 


.04-0.01 




A2199 


0, 


.20- 


-0, 


.30 


4, 


.18±0, 


.46 


0, 


.50±0, 


.13-0.01 




A496 


0, 


.00- 


-0, 


.03 


3 


.01±0, 


.02 


0, 


.93±0, 


.05+0.08 




A496 


0, 


.03- 


-0, 


.06 


4, 


.16±0, 


.03 


0, 


.70±0, 


.03+0.05 




A496 


0, 


.06- 


-0, 


.10 


4, 


.72±0, 


.05 


0, 


.61±0, 


.04+0.01 




A496 


0, 


.10- 


-0, 


.20 


4, 


.73±0, 


.08 


0, 


.48±0, 


.04+0.01 




A496 


0, 


.20- 


-0, 


.30 


3 


.94±0, 


.17 


0, 


.34±0, 


.05+0.03 




A1644 


0, 


.03- 


-0, 


.06 


4, 


.47±0, 


.43 


0, 


.76±0, 


.22+0.04 




A1644 


0, 


.06- 


-0, 


.10 


4, 


.89±0, 


.35 


0, 


.84±0, 


.17+0.01 




A1644 


0, 


.10- 


-0, 


.20 


4, 


.95±0, 


.30 


0, 


.64±0, 


.11-0.00 
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K. Matsushita: Fe abundance in the Intracluster Medium 



cluster r (riso) fcT" (keV) Fe (solar) from FfcHok/-F3.5_6 Fe (solar) from -FfoHk/-F3.5-6 
Best-fit^ ± A.^tat Af,, Best-fit^ ± A.^^.t 



A 1 f\dd 

rt.X 1^4:4: 


n 






u. 


.ou 


A 90+0 
4:. zunzu 


■^i 

. 1 




u. 


fil +0 
.u X mu. 


1 '=1+0 04 
. xo r U.U4 










n 


nn 




u. 


.uo 


A fiO-l-O 

4:.UZZIZU 


94 

. Z4: 


1 


17-1-0 
.o ( mu. 


1 7 01 

.J- ( U .U 1 










n 

u 


.yjo- 




-u. 


.uu 


A 71 -1-0 
^. ( imu 


90 
.zu 




u. 


QO-l-0 
.yumu. 


00+0 01 
.uy^^u-ui 










n 






u. 


. ±u 




90 
. zu 




u. 


. ( omu. 


00+0 09 
.uy n^u.uz 










n 


1 




u. 


.zu 


A fi'?-|-0 
4:.uoznu 


1 S 
. 1 o 




u. 


fiH-O 
.uomu. 


Ofi+0 09 
.uun^u.uz 








A '^'^(^9 


n 

u 


90 




-u. 


.ou 


A zlQ-l-0 
zcu 


.o^ 




u. 


.ozmu. 


14+0 09 

. ±4:T^U.UZ 








i\oooo 


n 


nn 

. uu 




u. 


.uo 


07-1-0 
o .u ( znu 


1 7 




u. 


. ( omu. 


07+0 07 

.U ( ^U.U ( 








i\oooo 


n 


.uo 




u. 


Ofi 
.uu 


o .ooznu 


1 7 




u. 


fi4-(-0 
.U4mu. 


07 01 

.U ( U .U 1 



u 


S1+0 
.oomu . 


1 Q 01 
. 1 y u.u X 


A '^'^'^R 


n 

u 


.uu- 




-u. 




o.o ( mu 


1 fi 

. lU 


u. 


^0+0 
.oumu. 


Ofi 01 
.uu-u.ui 



u 


40+0 
.4y mu . 


110 01 
, lO-U.Ul 


A'^'^'^S 
i\oooo 


n 


. lU 




u. 


90 
.zu 


o . z yznu 


1 1 




u. 


4fi-l-0 
.4umu. 


04+0 09 

.U4^U.UZ 



u 


^1+0 
. oomu . 


1 4 07 

. 14 U.U / 




n 


90 
. zu 




u. 


.ou 


A 77-1-0 
^.11 znu 


1 fi 
. lU 




u. 


41^ -1-0 
.4omu. 


0'=l+0 07 

.UOn^U.U ( 








A '^(^97 


n 

u 


00 
.uu- 




-u. 


.uo 


K 7fi-|-0 
o. ( umu 


"iO 
.oy 


n 

u. 


.oimu. 


00 00 
.uy-u.uu 








A'^B97 


n 


0"? 
.uo 




u. 


OR 
.uu 


o . 1 ^znu 


9Q 
. z y 




u. 


4Q-I-0 
.4y mu. 


07+0 01 
.u ( r u.u X 








A'^(^97 


n 


Ofi 
. uu 




u. 


1 n 

. ±u 


fi 1 1 -1-0 
u . 1 1 znu 


90 
. zu 




u. 


44+0 
. 44mu. 


Ofi 00 
.uu u .uu 








A '^(^97 


n 

u 


1 
. lu- 




-u. 


90 
.zu 


A Q7-I-0 
^.y ( mu 


1 7 




u. 


'=,9+0 
.ozmu. 


Ofi+0 Ofi 
.UUn^U.UU 








A'^fi97 


n 


90 
. zu 




u. 


.ou 


A 04-1-0 


44 

.4:4: 




u. 


71+0 
. / omu. 


1 7-1-0 1 7 








A 1 7Q^ 


n 


00 
. uu 




u. 


.uo 


A 91 -1-0 

4:. Z ' 1 U 


0^ 
.UO 




u. 


Sfi+0 
.oumu. 


04+0 1 1 

.U4^U. J- -L 








A 1 70"^ 


n 

u 


0'^ 




-u. 


Ofi 
.uu 


o.ozinu 


1 
. lU 




u. 


71+0 
. ( omu. 


Ofi+0 09 
.UUn^U.UZ 



u 


04+0 
. y4mu . 


91 04 
.ZO-U.U4: 




n 


Ofi 
. uu 




u. 


1 n 

. ±u 


o . yoznu 


1 4 




u. 


^fi+0 
.oumu. 


07+0 09 

.U ( T^U.UZ 



u 


^4+0 
. 04mu . 


1 4 Ofi 

, 14 U.UU 


A 1 7Q^ 


n 


1 

. lU 




u. 


90 
.zu 


fi 91 -1-0 

U . Z ' 1 u 


1 fi 

. lU 




u. 


41+0 
.4omu. 


0'=, 01 
.uo U .U 1 



u 


14+0 
. 04mu . 


OQ+0 00 
.uy T^u .uu 


A 1 70"^ 


n 

u 


90 
. zu- 




-u. 


"^O 
.ou 


ft 90-1-0 
u.zumu 


97 

.Z t 


u. 


'=,7+0 
.o ( mu. 


00+0 00 
.uy ^u.uu 








A 1 70"^ 


n 


"^0 
. ou 




u. 


•^0 
.ou 


K 79-1-0 
o . 1 zznu 


fi9 
. uz 




u. 


1'=l+0 
.oomu. 


1 0+0 01 
. ±un^u.uo 








iT-OrJ 


n 


00 
. uu 




u. 


0"^ 
.uo 


A 94-1-0 

4:. Z4:ZIZU 


1 ^ 
. 1 o 




u. 


QS+0 
. y omu. 


OQ+0 OQ 
.uy 1 u.uy 










n 

u 


0'^ 




-u. 


Ofi 
.uu 


K fiO-l-O 

o.uymu 


9'^ 
.zo 




u. 


84+0 
.o4:mu. 


1 1 01 










n 


Ofi 
. uu 




u. 


1 

. -LU 


o . c/zznu 


■^fi 
. ou 




u. 


fil +0 

.u -L mu. 


1 04 

. ±U U .U4 










n 


1 

. lU 




u. 


90 
.zu 


fi 1 1 -1-0 
u . 1 1 znu 


'^9 
. oz 




u. 


'=,9+0 
.ozmu. 


08 09 
.uo u .uz 










n 

u 


90 
. zu- 




-u. 


"^O 
.ou 


c; 97-1-0 
o.z ( mu 


^9 
.oz 


u. 


4Q+0 
.4y mu. 


1 9+0 0'=1 
. iz^u.uo 








iT-OUU ( 


n 


00 
. uu 




u. 


O"? 
.uo 


fi fil -1-0 
u .u 1 mu 


10 
.ou 




u. 


fil +0 
.u ± mu. 


1 00 
. ±u u .uu 




u 


S^+0 
. oomu . 


1 Q nil 
. 1 y U.UO 


iT-OUU ( 


n 


0'^ 
. uo 




u. 


Ofi 
.uu 


fi 41 -1-0 
u .4: ± mu 


1 fi 

. lU 




u. 


4Q+0 
.4y mu. 


07+0 01 

.U ( 1 u.u ± 



u 


^Q+0 
. o y mu . 


1 4 04 

. 14 U.U4 




n 

u 


Ofi 
.uu- 




-u. 


1 
. lU 


u.oomu 


1 4 


n 

u. 


'=11 +0 
.oimu. 


Ofi+0 09 
.uu^u.uz 



u 


fi7+0 
.u ( mu. 


1 1 0'=1 
. lO-U.UO 




n 


1 

. lU 




u. 


90 
.zu 


G Q9-I-0 
u . ozmu 


1 

. lU 




u. 


4S+0 
.4omu. 


Ofi+0 01 
.uu r u.u X 



u 


^1+0 
. oomu . 


1 01 

. lU U.UO 




n 


90 
. zu 




u. 


^0 
.ou 


u.oomu 


9Q 
. z y 




u. 


1fi+0 
.oumu. 


Ofi 04 

.uu U .U4 



u 


1Q+0 
.oy mu . 


OQ+0 09 
.uy T^u .uz 




n 

u 


■^o 

. ou- 




-u. 


.ou 


U.oomu 


9fi 
.zu 


u. 


9Q+0 
.zy mu. 


0'=,+0 01 
.UOn^U.Ul 








A49B 


n 


00 
. uu 




u. 


0"^ 
.uo 


Q 70-1-0 
o . 1 y mu 


01 

.U 1 




u. 


78+0 
. ( omu. 


09+0 1 1 

.UZ^U. -LO 










n 


0'^ 
. uo 




u. 


Ofi 
.uu 


o .oy mu 


04 

.U4: 




u. 


fifi+0 
.uumu. 


0'=l+0 14 

.UOn^U. -L4 



u 


70+0 
. 1 umu . 


1 ^ 90 
. 1 o u.zu 


A A9(^ 


n 

u 


Ofi 
.uu- 




-u. 


1 
. lU 


a 98-1-0 
u.zomu 


07 
.U ( 




u. 


'=14+0 
.04:mu. 


Ofi+0 Ofi 
.uu^u.uu 



u 


^fi+0 
.oumu. 


1 Ofi 
. lU-U.UU 




n 


1 

. lU 




u. 


90 
.zu 


u . 1 omu 


07 

.U 1 




u. 


^9+0 
.ozmu. 


07 01 

.U ( U .U 1 



u 


40+0 
.4umu . 


0^+0 01 

.uo T^U .U 1 




n 


90 
. zu 




u. 


^0 
.ou 


u . o^mu 


1 fi 

. lU 




u. 


48+0 
.4omu. 


07+0 01 

.U ( T^U.UO 



u 


11+0 
. oomu . 


OS 09 
.uo u.uz 


A49f^ 


n 

u 


■^o 

. ou- 




-u. 


.ou 


o.oomu 


94 

. Z4 




u. 


48+0 
.4:omu. 


Ofi+0 07 
.UUn^U.U ( 








A99'^fi 


n 


0'^ 
. uo 




u. 


Ofi 
.uu 


7 1 9-1-0 
( . 1 zmu 


fi'^ 
.UO 




u. 


'=l'=l+0 
.oomu. 


1 '=1 Ofi 

. -LO U .uu 








A99'^fi 


n 


Ofi 
. uu 




u. 


1 
. ±u 


a 70-1-0 
u . 1 umu 


.oo 




u. 


'=19+0 
.ozmu. 


1 0+0 09 
. ±u ru.uz 




u 


Sfi+0 
. oumu . 


91 a'\ 

. Z 1 U.UO 


A99'^(^ 


n 

u 


1 
. lu- 




-u. 


90 
.zu 


a 77-1-0 
u. ( ( mu 


.ou 


u. 


4fi+0 
.4:Umu. 


08 01 
.uo-u.ui 



u 


^9+0 
.ozmu. 


1 9 m 

. 1 Z-U.U-L 


A99'^fi 


n 


90 
. zu 




u. 


^0 
.ou 


o OS-l-0 
o .uomu 


41 




u. 


'=14+0 
.04mu. 


110 1^ 

. ±o U . 1 o 








A99'^fi 


n 


■^o 

. ou 




u. 


•^0 
.ou 


Q 1 1 -1-1 


9^ 
. zo 




u. 


41 +0 

.4± mu. 


1 fi 1 ^ 

. XU U . 1 o 








A '^'^71 


n 

u 


00 
.uu- 




-u. 


0"? 
.uo 


ft fifi-l-0 
u.uumu 


91 

. Z 1 




u. 


fi7+0 
.u ( mu. 


1 01 
.lU-U.Ul 



u 


fi^+0 
.uomu. 


1 1 09 
, lO-U.UZ 




n 


0'^ 
. uo 




u. 


Ofi 
.uu 


u . oomu 


91 
. z ± 




u. 


fi8+0 
.uomu. 


1 0+0 01 

. ±U T^U.U ± 



u 


84+0 
. 04mu . 


1 4 Ofi 

. 14 U.UU 


A '^'^71 


n 

u 


Ofi 
.uu- 




-u. 


1 
. lU 


ft 00-1-0 
u.yumu 


1 
. ly 




u. 


'=10+0 
.oy mu. 


00 0^ 
.uy-u.uo 



u 


^4+0 
.04mu. 


00+0 01 
.uy^^u.uo 


A '^'^71 


n 

u 


1 

. lU- 




-u. 


90 
.zu 


ft 71 -Lo 
u. ( imu 


9fi 
.zu 




u. 


'=lfi+0 
.oumu. 


00 09 
.uy-u.uz 



u 


10+0 
.oymu. 


00 01 
.uy-u.ui 


A'^^71 


n 


90 
. zu 




u. 


^0 
.ou 


7 89-1-0 
( .ozmu 


.oo 




u. 


fi1+0 
.uomu. 


910 19 

.Z ± U . -L Z 








A 9090 


n 

u 


00 
.uu- 




-u. 


0"? 
.uo 


ft zlS-l-0 
u.oomu 


1 7 


1 


08+0 
.uomu. 


1 fi+0 1 
. lUn^u. ly 


1 


91 +0 
. z imu. 


9fi 08 
.zu-u.uo 


A 9090 


n 

u 


0"? 
.uo- 




-u. 


Ofi 
.uu 


7 fiQ-i-O 
( .uomu 


90 
.zy 


u. 


8fi+0 
.oumu. 


1 7+0 09 
. 1 ( n^U.UZ 



u 


01+0 
.yomu. 


1 fi 01 

. lU-U.UO 


A 9090 


n 


Ofi 
. uu 




u. 


1 
. ±u 


7 zt ^-1-0 


.oo 




u. 


fi'=l+0 
.uomu. 


1 1+0 0'=1 

. ±0 T^U.UO 



u 


^7+0 
. o 1 mu . 


1 1 09 
. lo u.uz 


A 9090 


n 

u 


1 
. lu- 




-u. 


90 
.zu 


( . ( omu 


41 

.41. 


n 

u. 


'=17+0 
.o / mu. 


1 1+0 07 
. lO^U.U ( 



u 


7^+0 
. ( omu. 


1 4 01 
. 14-U.UO 


A 9090 


n 
u 


90 
. ZU- 


-u. 


"^O 

.ou 


o.uomu 


04 
. y4 


n 

u. 


44+0 

.4:4:mU. 


17 0^ 
. i ( -U.UO 








A1656 


0, 


.00- 


-0, 


.03 


8.17±0, 


.09 


0, 


.41±0, 


.07-f0.05 


0, 


.46±0. 


,03-0.02 


A1656 


0, 


.03- 


-0, 


.06 


8.22±0, 


.06 


0, 


.39±0, 


.07-0.03 





.41±0. 


,03-^0.01 
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.79±0. 


,17-0.04 
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cluster 


r (riso) 


kT"" (keV) 


Fe (solar) from FfcHok/-F3.5-6 
Best-fit^ ± A-,^,Ai, 


Fe (solar) from ^foHk/-F3.5-6 
Best-fit^ ± A,-t,t Af^3 








A3266 


0.03-0.06 


8.29±0.39 


0.58±0.13-0.02 


0.55±0. 11-0.00 


A3266 


0.06-0.10 


8.21±0.31 


0.49±0.10+0.00 


0.38±0.09-0.01 


A3266 


0.10-0.20 


8.09±0.21 


0.46±0.09-0.02 


0.46±0.06+0.01 


A3266 


0.20-0.30 


6.86±0.30 


0.33±0.07-l-0.03 


0.42±0. 12-0.00 


A3266 


0.30-0.50 


6.29±0.47 


0.48±0.12-h0.13 




A754 


0.03-0.06 


8.93±0.60 


0.45±0.12-h0.03 


0.62±0. 13-0.01 


A754 


0.06-0.10 


8.85±0.42 


0.58±0.13-0.02 


0.36±0.08-H0.01 


A754 


0.10-0.20 


9.66±0.38 


0.52±0.11-0.06 


0.47±0.06-f0.01 


A754 


0.20-0.30 


10.04±1.67 


0.38±0.16-0.04 


0.54±0. 12-0.01 



"Best-fit ICM temperature with la statistical error derived from the IT model fits 
''Best-fit Fe abundance using the ICM temperature derived from the IT model fits 

°1 a statistical error, considering a 10 % systematic uncertainty in the ICM temperature derived from the IT model. 
Difference in the derived Fe abundances using the best-fit multi-T model and the IT model 



